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ABSTRACT 

We report the discovery of blue shifted (Av>200km s _1 ) mid-infrared [Nelll] 
and/or [Nev] emission in 25 out of 82 ULIRGs (30% of our sample). The incidence 
of blue shifted [Nev] emission is even higher (59%) among the sources with a [Nev] 
detection — the tell-tale signature of an active galactic nucleus (AGN). Sixteen 
ULIRGs in our sample, eleven of which are optically classified as AGN, have 
[Nelll] blue shifts above 200 km s _1 . A comparison of the line profiles of their 
12.81 /xm [Nell], 15.56 //m [Nelll] and 14.32 /im [Nev] lines reveals the ionization 
of the blue shifted gas to increase with blue shift, implying decelerating outflows 
in a stratified medium, photo-ionized by the AGN. The strong correlation of the 
line width of the [Nelll] line with the radio luminosity indicates that interaction 
of expanding radio jets with the dense ISM surrounding the AGN may explain 
the observed neon line kinematics for the strongest radio sources in this sample. 

Subject headings: ISM: jets and outflows — infrared: ISM — galaxies: ISM - 
galaxies: active 



1. Introduction 

Galactic-scale outflows are now generally recognized to be important phenomena affect- 
ing or even regulating the evolution of galaxies. Starburst-driven winds, powered by the com- 
bined kinetic energy released in the outflows from massive stars and sup ernovae, are capable 



of depositing metals, dust and energy into the galaxy halo and beyond (jVeilleux et al.l 12005 
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and references therein). In galaxies hosting an active galactic nucleus (AGN), radiation- 
driven AGN winds and collimated outflows of radio plasma, or "jets" , offer another effective 
mechanism of transporting large amounts of energy and momentum out of the nucleus. In re- 
cent years, both proc esses have been impl i ed in disrupting/ fragmenting the obscuring cocoon 



in merger remnants (IHopkins et al.l 120051 ; iHolt et al. 



2008|). 



In ultraluminous infrared galaxies (ULIRGs; Lir=10 12-13 L ), molecular gas driven into 
the merging nuclei by the galaxy interaction is capable of supporting both a prolonged mas- 
sive starburst and/or feeding a (nascent) AGN. This makes ULIRGs prime candidates for 



observing extreme AGN and/or st arburst driven outflows. Indeed, observations o 



stellar Na I-D absorption line (e.g. lHeckman et al.ll200d ; lMartinll20 05; Ru pke et al. 



optical forbidden low a nd medium ioniza tion lines (e.g. iHeckman et al 



1990 



the inter- 



2005) and 



Wilman et al. 



19991 ; IHolt et al.l 120031 ; iLfpari et al.l 120031 ) indicate that outflows are ubiquitous in ULIRGs 



and that the mom e ntum and energy in jection is in most cases dominated by the starburst 
faupke et al.ll2005l : IVeilleux et al1l2005h . 

Mid-infrared kinematic studies of forbid den line emission i n AGN and ULIRGs are 
scarce. In their ISO-SWS survey of local AGN. ISturm et al.l (120021 ) reported the detection of 
line asymmetries and line shifts for some of their sources. Detailed studies could, however, 
only be performe d for the two brightest AGNs, NGC4151 Jsturm et alill999h and NGC 1068 
(ILutz et al.ll2000l ) . For these it was determined that dust entrained in outflows cannot explain 
the observed asymmetries in optical forbidden lines. With the advent of Spitzer-IRS sensitive 
medium- resolution spectroscopy has become available fo r kinematic studies of AGN and 
ULIRGs at a resolving power i?=600. iDasyra et al.l (120081 ) have exploited these capabilities 
and found the line widths of the high-ionization AGN narrow line tracers 14.32 /im [Nev] 
and 25.89 /im [Oiv] to be correlated with the black hole m ass for AGNs for wh ich the narrow 
line region is dominated by virial motion. More recently, ISpoon et al.l (120091 ) presented the 
discovery of strongly blueshifted ionized Ne + , Ne 2+ and Ne 4+ gas in three ULIRGs and 
proposed that the outflowing gas may be tracing the disruption of the obscuring medium 
around buried AGNs in these mergers. In this Letter we report the mid-infrared discovery 
of strongly blue shifted neon line emission in another 22 ULIRGs. 



Observations and data reduction 



The results prese nted in this paper a re bas ed on high resolut ion (i?~600; Av ~500km 
- 1 ) 10-37 /im Spitzer JWerner et alil2004h IRS Jrlouck et alil2004h spectra o f ULIRGs, AGN 



and starburst galaxies. The core sample is formed b y the IRS GTO sample (jWeedman et al. 



20051 ; iFarrah et al.l 120071 : iBernard-Salas et al.l 120091 ) which is not complete in any sense 
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All sources were observed in staring mode. The 10-19.5 /im portion of the spectra were 
obtained in the IRS Short High (SH) module and the 19.3-37/im sections in t he Long High 



(LH) module. The data reduction proceeded along the same steps as detailed in lSpoon et al. 



(120091 ). The wavelength calibration of our spectra should be accurate to about 1/5 (SSC 
priv. comm.) of a resolution element (500 km s _1 ), which amounts to 100 km s _1 . The 
accuracy is however known to vary from order to order and is generally lowest near the order 
edges and best toward the order center. 



3. Analysis 



We have used the non-linear least squares fitting IDL package MPFIT (jMarkwardtll2008l ) 
to perform single and double component gauss fits to the high-resolution Spitzer emission 
line spectra of the galaxies in our sample. With these fits we computed the full width and 
velocity center at 20% of the peak flux (FW20 and VC20, respectively) for all detected lines. 
Our choice of evaluating at the 20% level rather than at the 50% level is motivated by the 
presence of weak blue wings in the neon lines. 

In this particular study we focus on the mid-infrared fine structure lines of neon, because 
of their large range in ionization potential (IP = 21, 41, 97 and 127 eV for Ne + , Ne 2+ , Ne 4+ 
and Ne 5+ , respectively) and their superior detectability compared to fine-structure lines of 
argon and sulphur. Also, the rest wavelengths of the five neon lines are well away from the 
strong Si-0 stretching resonance of amorphous silicate grains at 8-12/im which is responsible 
for a factor 3-5 increase in opacity compared to surrounding wavelengths. Furthermore, the 
three strongest lines, 12.81 /im [Ne n], 15.56 /im [Ne ill] and 14.32 /xm [Ne v], all fall within 
the same IRS slit, SH, for redshifts up to 0.25. At higher redshifts the 15.56 /im [Ne ill] line 
shifts from the SH slit into the LH slit, while above z=0.32 the 7.65 /im [Ne vi] line shifts 
into the SH slit. Our neon line sample consists of 82 ULIRGs, 3 HyLIRGs, 33 AGNs (both 
Seyferts and QSOs) and 15 starburst galaxies suited for detailed line profile analysis. 

In Fig.[T]we plot for the [Ne ill] and [Ne v] lines the FW20 as a function of VC20 for 
all sources in our sample. The panels strikingly illustrate the predominance of blue over 
red shifts among the neon line detections, and the tendency for the VC20 blue shift to be 
correlated with the FW20 line width. The latter appears to be best established for the [Ne 
ill] line. 

Among our ULIRG sample 25/82 sources (30%) show blue shifted [Ne ill] and/or [Ne 
v] emission. This fraction excludes sources with VC20 blue shifts less than 200 km s _1 , 
for which we cannot rule out uncertainties in the order-to-order wavelength calibration to 
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mimick a real blue shift. 21/25 ULIRGs (84%) are optically classified as AGN. Among the 
32/82 ULIRGs with a [Ne v] detection, the [Ne v] line is found to be blue shifted (VC20<- 
200km s _1 ) for 19 ULIRGs (59%), 6 of which are optically classified as type- 1 and 10 as 
type-n AGN. 

In Fig. [T] the broadest [Nev] lines not showing clear indications for blue shifts have FW20 
line widths of 1500 km s -1 . Assuming these l ines to arise i n NLR gas in keplerian rotation 



around the central black hole, equation 1 of iDasyra et al.l (120081 ) predicts their black hole 

I - 



masses to be of the order of 10 9 ' 3 Mm. 



In Figs.[2]and[3]we show overlays of the velocity profiles of the mid-infrared neon lines for 
16 sources which have [Nelll] VC20 line center blue shifts of 200 km s _1 or more. The sources 
are sorted by decreasing FW20 of the [Nelll] line. Eleven of these are optically classified as 
AGN. The 16 panels reveal a great diversity in line shapes and line center shifts, both from 
line to line for a single target, as well as from object to object. In some sources the VC20 
line center shift is the result of a pronounced blue wing to the systemic component, while 
in others also the line peak is blue shifted. The latter is more often the case for the [Ne v] 
line than for the [Ne ill] line, and can also be more extreme. For other sources two gaussian 
components can be discerned, one at systemic velocity, the other at blue shifted velocities. 
Examples are the [Ne n] line of IRAS 07598+6508 and the [Ne ill] line of IRAS 15462-0450. 
In Tabled] we list the measurements of the kinematic parameters for the sources in Figs. [2] and 
[3] along with their computed 8-1000 /im infrared luminosities and 1.4 GHz radio luminosities. 

With lines of various ionization stages of neon measured in the same aperture, it is 
possible to investigate how the ionization of the gas changes as a function of velocity relative 
to systemic. In Figs. [2] and [3] we show the change in [Ne lll]/[Ne n] and [Ne v]/[Ne ill] line 
ratios at velocity intervals of 750 km s _1 with respect to systemic. Results are only shown 
for those velocity bins for which the ratio-ed lines are at least 1/20 of the peak line flux 
throughout the central 600 km s _1 of the velocity bin. For nearly all the sources for which 
either line ratio is defined over more than one velocity bin, the ionization appears to increase 
with increasing blue shift. In some cases this can be substantial, like for IZwl, where the [Ne 
v]/[Ne ill] ratio increases from 0.21 at systemic velocity to 1.3 in the velocity bin centered 
at -1500km s -1 . A clear exception is IRAS F00183- 7111, for which the [Ne lll]/[Ne il] ratio 
is remarkably constant over a wide velocity range JSpoon et al.ll2009h . The above trends 



can also be recognized in a more visual way by comparing for each source which neon lines 
have the most pronounced blue wings. This reveals that for the large majority of sources in 
Figs.[2]and[3]the blue wing of the [Ne v] line extends beyond that of the [Ne ill] line, while 
the blue wing of the [Ne ill] line extends beyond that of the [Ne n] line; consistent with an 
increase in ionization with increasing blue shift. 
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All sources in Figs. [2] and [3] display an unresolved or barely resolved [Ne n] component 
at systemic velocity. Some of these sources also have a [Ne ill] or even [Ne v] systemic com- 
ponent. The excitation of this neon gas holds clues as to the nature of the dominant nuclear 
ionizing source. In pure starburst galaxies the excitation of the nuclear starbu rst as measured 



from the [Ne Hi] / [Ne ll| line rati o is usually low and in the range of 0.05-0.2 (IThornley et al. 



200(1 iBernard-Salas et all liooiM In our sample of sixteen excitation conditions like these 



are met for the systemic components of IRAS 05024-1941, IRAS 07598+6508, IRAS 12127- 
1412NE, IRAS 15462-0450 and IRAS 15225+2350. Among these IRAS 15225+2350 is the 
only source optically classified as a starburst. For other sources, like IRAS 12514+1027 and 
IRAS 16156+0146NW, the [Ne m]/[Ne il] ratios for the velocity component at systemic ve- 
locity are far higher, 0.97 and 0.78, respectively. In these galaxies the AGN contributes 
to the systemic [Ne ill] emission, given the detection of a systemic [Ne v] emission compo- 
nent. For IRAS 05189-2524 and IRAS 11119+3257 the clear blue shift of both the [Nem] and 
[Nev] emission implies that none of the higher ionization gas in these sources has a keplerian 
component. 

In Fig.H] we plot the FW20 line width of the [Ne ill] and [Ne v] lines as a function 
of 1.4 GHz radio luminosity for all AGNs and ULIRGs in our sample. This reveals a clear 
trend of increasing line width with increasing r a dio po wer. A simila r trend was first found 
for the 5007 A [O ill] line by Iwilson fc Wfilisl Jl98ph and Iwhittle! Jl985h for samples of 
radio quiet (L[1.4GHz]<10 23 W Hz -1 ) seyfert galaxies and QSOs. In our data the trend 
appears to persist - albeit with some outliers - well into the radio-loud domain. The sources 
with strong blue shifts in their [Nelll] lines, all labeled in red, appear to define the upper 
envelope of this trend or a separate one. This is especially clear in the upper panel. The 
only obvious exception is the prototypical narrow-line Seyfert 1 (NLS1) galaxy IZwl. Note 
that IRAS 07598+6508 is not included, as it lacks a 1.4 GHz flux measurement. 



4. Discussion 

A striking result of our study is the identification of a sizeable subsample of sources 
which show a strong blue shift of the [Nev] line center in combination with a smaller blue 
shift of the [Nelll] line center (Figs. [2] and [3]) in their spectra. A per-source comparison 
of the [Nell], [Nelll] and [Nev] profiles strongly suggests these shifts to be the result of 
an increase in ionization with increasing blue shift in the line. This is supported by the 



1 A notabl e exception may be the pu re starburst galaxy NGC 7714, which has a measured [Ne m]/[Ne n] 
ratio of 0.75 (jlBernard-Salas et alJliooj ) 
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behavior of the [Nein]/[Neii] and [Nev]/[Nem] line ratios as a function of line velocity. The 
most extreme sources in this subsample do not even exhibit a [Nev] component at systemic 
velocity Examples are IZwl, IRAS 05189-2524, IRAS 07598+6508, IRAS 11119+3257 and 
IRAS 15462-0450. This strongly suggests that at least in these sources the observed [Nev] 
emission cannot arise in NLR gas in keplerian rotation. Instead, the observed line shifts 
and blue wings imply a nucle ar outflow in a hig hly stratified ISM, ionized by the central 
source (e.g. IRAS 13451 + 1232; ISpoon et al.ll2009i ). The observed increase of ionization with 
increasing blue shift then indicates that the bulk of the observed neon gas is moving toward 
us, and that the outflow speed decreases with distance to the ionizing source. The absence 
of a matching redward line component to the [Nev] line implies that the line of sight to the 
far side of the outflow is blocked at mid-infrared wavelengths. 

Note that most of these sources are also known from optical spec t roscop y to exhibit 
outflows. Nine of the 16 sources in Tabled] are listed by iLipari et al.l (120031 ) as showing 
optical outflow signatures, but also IZ wl and IRAS 23060+0505 have known optical outflows 
jLaor et al.lll997l : Iwilman et al.lll999h . 



Res ults similar to ours were recently reported for a sample of "blue outlier" NLS1 
galaxies ( Komossa et al. For these galaxies a similar correlation seems to exist between 

blue shift and line width of the 5007 A [Olll] line as we see for the 15.56 /iin [Nem] line (and 
for the 14.32 /im [Nev] line to a lesser degree). The authors favor a scenario in which gas 
clouds from the narrow line region are entrained in a decelerating wind, powered either 
by collimated radio plasma jets, or radiation pressure resulting from accretion close to or 
above the Eddington limit. Similar [Olll] blue shifts and line w idths have also been rep orted 
for two young dust -enshrouded rad io galaxies: PKS 1549-79 (ITadhunter et al. 1 120011 ) and 
IRAS 13451+1232 (IHolt et al.ll2003l ). Both are late stage mergers and the latter is one is 
part of our sample (Tabled]). The broad, blue shifted [Olll] emission from these sources 
may trace the ex pansion of radio jets through the obscuring shell around the buried AGN 
(|Holt et allboosh . 



Given the high radio luminosities for IRAS 13451+1232, IRAS F00183-7111, IRAS 11119+3257, 
IRAS12127-1412NE, and IRAS 05024-1941, the apparent correlation between radio lumi- 
nosity and [Nein] line width in Fig. [4] may indicate that for these sources interaction of 
expanding radio jets with the dense ISM surrounding the AGN may explain the observed 
neon line kinematics. 

Four of the sources in Fig. [2] which show outflow signatures in their neon lines do not 
show the tell-tale mid-infrared signature of AGN activity, the presence of a 14.32 /iin [Nev] 
line, in their spectrum. While this may indicate their outflows to be starburst-driven, three 
of these sources, IRAS F00183-7111, IRAS 05024-1941 and IRAS 12127-1412NE, display 



- 7- 



strong evidence for the presence of a powerful AGN - e.g. 1.4 GHz radio fluxes 4-10 times 
in excess of the value expected for star forming galax ies. In addition, I RAS F00183-7111 
and IRAS 12127-1412NE have silicate strengths (S sU ; ISpoon et all I2007T I below -2.5, indi- 
cating an apparent silicate optical depth of at least 2.5 toward the nucleus. Hence, it is 
conceivable that the non-detection of [Nev] emission in these sources is the result of ob- 
scuration (along our line of sight) of the interior, most highly ionized part of the outflow 
rather than the absence or insignificance of an AGN in these sources. For the fourth galaxy, 
IRAS 01003-2238, conflicting classi fications exist . Several studies have conclud e d that this 



galaxy may be starburst -powered (1 Armus et all 1 19881 ; iRigopoulou et all Il999l ; llutz et al. 



1999; 



Wilman et all 119990 . whereas lAUen et all (Il99ll ) classify this galaxy as Seyfert-II and 



lipari et a 



(20 03j) refer to it as a QSO. Its Spitzer-IRS spectrum is continuum-dominated 



(class 1A; ISpoon et all 120071 ) with a 6.2 /zm PAH equivalent width more than ten times 
lower than for starburst galaxies. From the PAH lu minosity we calcula te a contribution of 
the starburst to the infrared luminosity of 10-20% (IPeeters et all 12004 ). This would imply 
that IRAS 01003-2238 is AGN-dominatecfl and its outflow likely AGN-powered. The mere 
detection of blue wings in [Nell] and [Nein] line profiles in Spitzer-IRS-SH galaxy spectra 
may hence be a general signpost of AGN activity. 



The author thanks Ian Waters for help with measuring the line profiles and Jeremy 
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Table 1. Neon emission line kinematics 



Target name 
(IRAS) 


z a 


class b 


FW20 [NeII] 
(km s _1 ) 


FW20 [NeIII] 
(km s" 1 ) 


FW20 [NeV ] 
(km s -1 ) 


VC20[Ncii] 
(km s" 1 ) 


VC20 
(km 


[NcIII] 
s" 1 ) 


VC20[n o v] 
(km s- 1 ) 


log L(IR) C 
(Lq) 


log L(1.4GHz)' 
(W Hz- 1 ) 


00183-7111 


0.3219 


L/A (1,2) 


3037 


2807 


- 


-393 




-449 


- 


12.93 


26.04 (1) 


IZwl 


0.06089 


NLS1 (3) 


<400 


2543 


2524 







-988 


-1190 


11.95 


22.86 (2) 


01003-2238 


0.1179 


H/S2 (4,7) 


<400 


1837 


- 







-559 


- 


12.29 


23.68 (2) 


05024-1941 


0.1937 


S2 (5) 


773 


2026 


- 


-66 




-446 


- 


12.39 


24.49 (2) 


05189-2524 


0.04274 


S2 (5) 


580 


1189 


1741 







-366 


-1032 


12.16 


23.01 (3) 


07598+6508 


0.1484 


Sl c (5) 


2109 


2678 


1073 


-626 




-469 


-847 


12.53 




11119+3257 


0.1903 


SI (5) 


473 


2226 


1677 







-679 


-870 


12.66 


25.04 (2) 


12071-0444 


0.1286 


S2 (5) 


<400 


1155 


945 


-23 




-320 


-350 


12.39 


23.52 (4) 


12127-1412NE 


0.1330 


L (5) 


986 


1794 




-146 




-484 




12.17 


24.56 (2) 


12514+1027 


0.3194 


S2 (6) 


969 


977 


1756 







-203 


-387 


12.77 


24.44 (2) 


13451+1232 


0.1214 


S2 (5) 


1228 


2841 


3072 


-109 




-821 


-1340 


12.29 


26.31 (2) 


15130-1958 


0.1103 


S2 (5) 


702 


1721 


1281 







-496 


-535 


12.14 


23.46 (2) 


15225+2350 


0.1386 


H (5) 


<400 


1476 


1154 


-14 




-266 


-566 


12.16 


23.44 (4) 


15462-0450 


0.1001 


NLS1 (3) 


<400 


1703 


634 


-16 




-403 


-1090 


12.20 


23.44 (4) 


16156+0146NW 


0.1326 
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Note. — Results from the 1 or 2-component gauss fits to the 12.81 [Nell], 15.56 [Nein] and 14.32 [Nev] line profiles. FW20 is defined as the full width at 20% 
of the peak flux, and is corrected for the instrumental profile. We estimate the uncertainty in these quantities to be 5, 10 and 15% for the [Nell], [Nein] and [Nev] 
line, respectively. VC20 is the velocity shift of the line center measured at 20% of the peak flux. The uncertainty in the VC20 values of the neon lines are limited 
to order-to-order wavelength calibration uncertainties which we estimate to be smaller than 150 km s . 

a The adopted redshift (z) is defined by the line center of the [Nell] line measured at 80% of the peak flux (VC80). 

b Optical spectral classifications: (SI) S eyfer t type 1, (NLS1) narrow -line Seyfert type 1, (S2) Se y fert ty pe 2, (L) LINER, (H) Hll, (A) Ambiguous. Classifica tions 
were taken from (1) lArmus et alj lll989h, ( 21 iBuchanan et al.l ll2006ll , (3) IVeron-Cettv fe Veronl ll2006r i. (4) IVeilleux et all Jl999ll , (5) IVeilleux et ahl (l2002h . (6) 
IWilman et al] <1999l) . (7l lAllen et al.l lll99lll 

C (IR)=L (8-1000 /im) as defined bv lSanders fc Mirabell lll996ll 

d Radio luminosities were compiled from: (ll lDrake et~afl (l2004h . (2) NVSS llCondon et al.lll998ll . (3l lYun et ahl J200lll . (4) FIRST l lWhite et al.lll997ll 
e Broad absorption line (BAL) QSO l lLi'parilll994ll 
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Fig. 1.— The full width at 20% of the peak flux (FW20) for the line profiles of 15.56 fj,m [Ne 
ill] and 14.32 fim [Ne v] lines shown as a function of the shift of the line center measured 
at 20% of the peak height (VC20). Typical uncertainties for the FW20 values are 10% and 
15% for the [Ne ill] and [Ne v] lines, respectively. Upper limits on FW20 values denote 
unresolved lines. The sources are color-coded as follows: ULIRGs are red, starburst galaxies 
blue, QSOs and Seyfert-1 galaxies light green, and Seyfert-2 galaxies dark green. 
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Fig. 2. — Comparison of the line profiles of 12.81 /um [Ne n] (red), 15.56 /mi [Ne ill] (blue) 
and 14.32 /mi [Ne v] (green), as observed with the Spitzer-IRS high-resolution spectrographs 
at it!=600. The sources are sorted in order of decreasing FW20 of the [Nelll] line. Systemic 
velocity is indicated by vertical dashed lines. Data points are shown as diamonds, fits to the 
data points are shown as continuous lines. The filled gray area in the top-left panel shows 
the gauss profile of an unresolved emission line at a resolution of i?=600. Also plotted are 
the line ratios of [Ne lll]/[Ne il] (gray open circles) and [Ne v]/[Ne ill] (gray open squares) 
at 750 km s _1 velocity intervals. These ratios are computed from the fitted line profiles. The 
number in the upper right corner is the VC20 line shift of the [Ne ill] line in km s _1 . 
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Fig. 3. — Comparison of the line profiles of 12.81 /im [Nell] (red), 15.56 /mi [Nelll] (blue), 
14.32 //m [Nev] (green) and 7.65 /mi [NeVl] (brown; only IRAS 12514+1027), as observed 
with the Spitzer-IRS high-resolution spectrographs at i?=600. The lay-out is the same as in 
Fig.® 
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Fig. 4.— The full width at 20% of the peak flux (FW20) for the 15.56 fxm [Ne ill] and the 
14.32 //m [Ne v] line profile shown as a function of the 1.4 GHz radio luminosity density. 
Typical uncertainties for the FW20 values are 10% and 15%, respectively. Upper limits on 
FW20 values denote unresolved lines. The sources shown in red all have [Ne ill] VC20 blue 
shifts of 200 km s -1 or higher and are the same sources as listed in Tabled] 



